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; Abstract 

We use recent lattice data on the gluon and ghost propagators, as well as the Kugo-Ojima 
function, in order to extract the non-perturbative behavior of two particular definitions of the 

^: 

QCD effective charge, one based on the pinch technique construction, and one obtained from the 
standard ghost-gluon vertex. The construction relies crucially on the definition of two dimensionful 
quantities, which are invariant under the renormalization group, and are built out of very particular 
combinations of the aforementioned Green's functions. The main non-perturbative feature of both 
effective charges, encoded in the infrared finiteness of the gluon propagator and ghost dressing 
function used in their definition, is the freezing at a common finite (non-vanishing) value, in 
agreement with a plethora of theoretical and phenomenological expectations. We discuss the 
sizable discrepancy between the freezing values obtained from the present lattice analysis and the 
^ ■ corresponding estimates derived from several phenomenological studies, and attribute its origin to 

the difference in the gauges employed. A particular toy calculation suggests that the modifications 
induced to the non-perturbative gluon propagator by the gauge choice may indeed account for the 
observed deviation of the freezing values. 

PACS numbers: 12.38.Lg, 12.38.Aw 12.38. Gc 
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I. INTRODUCTION 



In recent years, a large number of independent lattice simulations have furnished highly 
non-trivial information on the infrared (IR) behavior of two fundamental ingredients of pure 
Yang-Mills theories, namely the (quenched) gluon and ghost propagators, for both SU(2) 
and SU (3) lH5|. In particular, these simulations have firmly established that (in the Landau 
gauge) the QCD gluon propagator and the ghost dressing function are IR finite and non- 
vanishing. Given that the entire issue is under intense scrutiny, it is natural to explore some 
of the most salient theoretical and phenomenological implications of these lattice results. 
The purpose of the present work is to use the available lattice data to extract the running 
of the QCD effective charge for a wide range of physical momenta, and, in particular, its 
behavior and value in the deep IR. This quantity lies at the interface between perturbative 
and non-perturbative effects in QCD, providing a continuous interpolation between two 
physically distinct regimes: the deep ultraviolet (UV), where perturbation theory is reliable, 
and the deep IR, where non-perturbative techniques must be employed. 

The generalization of the concept of the renormalization group (RG) invariant and process 
independent effective charge from QED to QCD is far from obvious, and has been discussed 
extensively in the literature. In this article we will consider two of the most standard defini- 
tions of the QCD effective charge. The first charge, to be denoted by a(q 2 ), constitutes the 
most direct non-Abelian generalization of the QED effective charge. This charge is obtained 
within the framework of the pinch technique (PT) 6 s| , and its generalization, known as 
generalized PT (GPT), introduced in [9|. Of particular importance in this construction is the 

nn 

profound correspondence [10|, lU| that exists between the PT (GPT) and the background- 
field method (BFM) [l^. The second charge, to be denoted by a^q 2 ), involves the ghost 
and gluon self-energies, in the Landau gauge, and in the kinematic configuration where the 
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14| becomes applicable |l5| . 



well-known Taylor non-renormalization theorem 

Both effective charges mentioned above display a strong dependence on the detailed 
characteristics of some of the most fundamental Green's functions of QCD. Specifically, in 
the case of a gh (q 2 ) the required ingredients are the conventional gluon propagator, A(g 2 ), 
(that of the gauges) and the ghost dressing function, F(q 2 ); both quantities are simulated 
on the lattice, and we will use them as inputs for obtaining a gh (q 2 ). For a(q 2 ) the situation is 
slightly more involved. The fundamental ingredient one needs for obtaining a(q 2 ) is the gluon 
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propagator of the PT-BFM, denoted by A(g 2 ), which, unfortunately, has not been simulated 
on the lattice yet. The way one establishes the required connection between the conventional 
gluon propagator (simulated on the lattice) and the PT-BFM propagator entering into the 
definition of a(q 2 ) is by resorting to two powerful non-perturbative identities. First, a formal 
relation known as "background-quantum" identity [ly, [l7|, given in Eq. (12.10p . relates the 
two gluon propagators by means of a special function, G(q 2 ), which plays a central role 
in the new Schwinger-Dyson equations (SDE) derived within the PT framework fis| . In 
fact, interestingly enough, in the Landau gauge only, G(q 2 ) coincides with the well-known 
Kugo-Ojima function 19]. The second identity, given in Eq. (I2.15p . allows one to obtain 
G(q 2 ) from F(q 2 ), to a very good approximation, given that the function that controls their 
difference, L(q 2 ), is numerically rather small, and vanishes exactly at q 2 = 0. Therefore, 
even though the theoretical origin of these two effective charges is vastly different [e.g., a(q 2 ) 
originates from a propagator, while a gh (g 2 ) from a vertex], they are very close in the entire 
range of physical momenta, and exactly coincide in the deep IR 20]. 

A large number of theoretical and phenomenological studies, based on a-priori very dis- 
tinct approaches (tI I21I 32 1 support the notion of the "freezing" of the QCD running coupling 
in the deep IR. In fact, when the QCD charge is constant (non- vanishing!) in the IR, and 
the quark masses are ignored, QCD becomes conformally invariant. Therefore, as has been 
emphasized amply in the recent literature [33J, the IR finiteness of the QCD effective charge 
constitutes a crucial requirement for the applicability of the powerful AdS/CFT correspon- 
dence [34j. 

As has been argued in numerous works, the IR finiteness of the effective charge and that of 
the gluon propagator are inextricably connected: they can be both traced back to the same 
phenomenon, namely the non-perturbative generation of a dynamical gluon mass, through 
the implementation of the Schwinger mechanism at the level of the SDE governing the gluon 
propagator 35]. Within the PT-BFM framework, the SDE solutions for the gluon self- 
energy, denoted by A(g 2 ), are used to form the RG-invariant combination d(q 2 ) = g 2 A(q 2 ) 
which, in turn, may be cast in the form d~ l (q 2 ) = [q 2 + m 2 (g 2 )]{feln( <? +4 ^ ^ q - 1 )}, where b 
is the first coefficient of the QCD /3 function, and A the QCD mass scale of a few hundred 
MeV. The non-perturbative generalization of a(q 2 ), the QCD effective charge, is contained 
in the curly brackets; evidently, the m 2 (q 2 ) in the argument of the logarithm tames the 
Landau pole, and a(q 2 ) freezes at a finite value in the IR, namely a -1 (0) = 61n(4m 2 (0)/A 2 ). 



3 



The IR finiteness of the effective charge obtained from the lattice data becomes manifest in 
the following way. First, one uses the available data for the gluon, the ghost, and the Kugo- 
Ojima function, to construct the lattice version of the corresponding dimensionful (mass 
dimension —2 ) RG- invariant quantity, denoted by d(q 2 ) in the case of a(q 2 ) (as above), 
and r(q 2 ) in the case of a gh (q 2 ) [see Eq. (I2.13P and Eq. (12. 3p . respectively]. The next step 
is to extract from d(q 2 ) and r(q 2 ) a dimensionless quantity, which will correspond to the 
associated effective charge. Both RG-invariant quantities have the the gluon propagator, 
A(g 2 ), as a common ingredient. Given that A(g 2 ) is effectively massive in the IR, one 
should follow the standard procedure used for massive gauge bosons, such as the W and 
the Z, namely factor out a massive "tree-level" propagator of the form [q 2 + m 2 (q 2 )]^ 1 . 
The procedure outlined above guarantees the freezing of the resulting coupling at a finite 
(non-vanishing) value. If, instead, a q~ 2 is factored out of the IR finite gluon propagator, 
one obtains (trivially) an effective charge that vanishes in the IR as q 2 . 

The article is organized as follows. In Section II we briefly review the definitions of 
the two effective charges under study, and recall the fundamental identities, Eq. ( 12 .101) 
and Eq. (12.151) . which relate their ingredients. Section III contains the main results of 
this work. In particular, after reviewing some of the most important lattice results on the 
(Landau gauge) gluon and ghost propagators, we construct the QCD effective charges and 
determine their freezing value in the deep IR. In section IV we discuss the sizable discrepancy 
between the freezing values obtained in the previous section and those favored by a variety 
of phenomenological studies. We argue that the main reason for the observed discrepancy 
is the difference in the gauge used: while the a(0) extracted from the lattice corresponds to 
the BFM Landau gauge, the phenomenological constraints are almost exclusively obtained 
in the BFM Feynman gauge. We will then derive an approximate formula that relates the 
two, suggesting that the discrepancy may be indeed accounted for by difference in gauge 
choices. Finally, in Section V we will discuss our results and present our conclusions. 

II. THE TWO EFFECTIVE CHARGES: DEFINITIONS AND BASIC CONCEPTS 

Before introducing the definitions of the effective charges and some of the important con- 
cepts related to them, we establish the necessary notation. The gluon and ghost propagator 
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will be defined as 



A M „(g) 



(2.1) 



D(g 2 ) = ^M, (2.2) 
g 2 

where £ denotes the gauge-fixing parameter, and P^ v {q) = g^v — q^qv/q 2 is the usual trans- 
verse projector. One has A _1 (g 2 ) = q 2 + zll(g 2 ), with II^(g) = P M1/ (g)II(g 2 ) the gluon 
self-energy; finally F(q 2 ) is the so called ghost dressing function. 

A reasonable definition of the QCD effective charge may be obtained from the ghost-gluon 



vertex in the Landau gauge [15l. |23|. Exploiting the fact that, in this gauge, the ghost-gluon 



vertex does not get renormalized, one can construct the RG-invariant product 

f(g 2 )=sV)A(g 2 )P 2 (g 2 ). (2-3) 

From this quantity one defines the effective charge as 

a gh (g 2 ) = [g 2 + m 2 (g 2 )]f(g 2 ), (2.4) 

where a(/i 2 ) = g 2 (/i 2 ) / '47r . Since A(g 2 ) and P(g 2 ) corresponds exactly to the quantities 
measured directly on the lattice, this definition constitutes the most direct way of extracting 
the non-perturbative QCD charge from the lattice. It should be noted, however, that away 
from the Landau gauge additional information on the form-factor of the ghost-gluon vertex 
must be supplemented, in order to define the RG-invariant quantity analogous to the r(g 2 ) 
of (I2.3p . This necessity, even though is not a limitation of principle, brings about several 
ambiguities; for example, the aforementioned vertex form-factor depends on two physical 
momenta, and a particular choice of the scale must be implemented, in order for the effective 
charge to be a function of a single momentum scale. In other words, one cannot obtain a 
universal definition of the charge, i.e., one that does not depend on the specific kinematic 
details of the vertex employed. 

A universal (process-independent) definition of an effective charge for every gauge may be 
obtained from the gluon self-energy in the (covariant) BFM. As is well-known, this quantity, 
to be denoted by fl^ Q )(g), captures the running of the QCD coupling for every value of the 
(quantum) gauge-fixing parameter, £ Q . In particular, at one loop, we have 

iU^(q) = q 2 g 2 [b\n (-g 2 /^ 2 ) + QJ, (2.5) 
5 



where b = 11Ca/487t 2 is the first coefficient of the QCD (3 function (j3 = —bg 3 ) in the 
absence of quarks, Ca is the Casimir eigenvalue of the adjoint representation (Ca = N for 
SU(N)), and the gauge-dependent constant C^ Q is given by (third item in (l0() 



Ca 



° iQ 16tt 2 



U(7 + U 67 



(2.6) 
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Note that the value £ Q = 1, i.e., the Feynman gauge of the BFM, is very special, because 
it reproduces the (gauge-independent) PT gluon self-energy; in this privileged gauge all 
unphysical longitudinal terms appearing inside an ostensibly gauge-independent quantity 
(physical on-shell amplitude, Wilson-loop, etc) have been discarded. 

For asymptotically large momenta one may neglect the constant C^ Q next to the leading 
logarithm, and write, in any gauge, (Euclidean momenta) 

^ (<?2) = q 2 [l + bg 2 ln(g 2 //i 2 )]' (2 ' 7) 
It is then easy to establish (e.g., by resorting to the QED-like identity = Z g , valid in 

the BFM to all orders and for every £ Q [12j) that the product 

d{q 2 )=g 2 ^ 2 )A{q 2 ), (2.8) 

is invariant under the renormalization group, i.e., it is an RG-invariant quantity, just as the 
r(q 2 ) in (12. 3p . From d(q 2 ) one may extract the QCD effective charge exactly as in (12. 4p . 
namely 

a(q 2 ) = [q 2 + m 2 (q 2 )]d(q 2 ). (2.9) 

In order to make contact between A(g 2 ) appearing in the definition of the RG-invariant 
product d(q 2 ) and the conventional propagator A(g 2 ) simulated on the lattice (in the Landau 
gauge), we employ a formal all-order identity, which relates them as follows [lfj Q 

A(q 2 )=[l + G(q 2 )] 2 A(q 2 ). (2.10) 

In the above formula the two gauge fixing constants, £ and £ Q , associated with A(q 2 ) and 
A(g 2 ), respectively, must be equal (but otherwise arbitrary); in particular, in the Landau 
gauge, £ = £ Q = 0. 

The function G(q 2 ) appearing in (12. lOj) is the g^ v component of a particular two-point 
function, denoted by A AtJ/ (g), defined as 

KM) = ~W 2 Ca [ H^D(k + q)A^(k)H av (k,q) 
Jk 

= ^G(g 2 ) + M^(g 2 ), (2.11) 
1 



where f, = /i e (27r)~ d J d d k, with d = 4 — e the space-time dimension. The function H„ v {k, q) 
appears in the all-order Slavnov- Taylor identity satisfied by the three gluon vertex 36j, and 
is related to the ghost-gluon vertex T^k, q) through the identity 



q v H lxv {k,q) = -iTJk,q). 



(2.12) 



At tree level, HjS) = ig^ v and T^(k,q) = T^(k,q) = —q^. Note that both G(q 2 ) and L(q 2 ) 
depend explicitly on the value of the gauge-fixing parameter £. 

Since the origin of the identity in (I2.10p is the BRST symmetry of the theory, it does not 
get deformed by the renormalization procedure. Thus, one can write (12. 8 j) in terms of A(g 2 ) 
and G(q 2 ) as follows 

Stf) = (2,3) 



-.(0) 



It turns out 



[1 + G(q 2 )] 



37H39J that the function G(q 2 ), coincides (in the Landau gauge only) with 



n 



the well-known Kugo-Ojima function [19(, u(q 2 ) defined as 



^S^ + P^q^uiq 2 ), 



(2.14) 



where (D^$) m = d^^ m + gf mnr A™® r is the usual QCD covariant derivative. The Kugo- 
Ojima function has been simulated on the lattice by means of Monte-Carlo averages of 
the operator time-ordered product appearing on the left-hand side of the defining equa- 
tion (I2.14p [del ]. Given that G(q 2 ) = u(q 2 ), the lattice information on u(q 2 ) may be used, in 
principle, into (12 . 1 3j) . together with the lattice results for the Landau gauge A(g 2 ). 

In addition, G(q 2 ) is related to the ghost dressing function F(q 2 ) and the form-factor 
L(q 2 ) of f l2TTTj) through the BRST identity 



l + G{q 2 ) + L{q 2 ). 



(2.15) 



This identity, in conjunction with the corresponding dynamical equations 20[ given in 
Eq. (13.51) . allows the indirect determination of G(q 2 ) and L(q 2 ) from the lattice data on 
the ghost dressing F(q 2 ) [41] . Thus, provided that one carries out the renormalization 
procedure in a way that manifestly preserves (12.151) . the two effective charges are related 
through the equation [20 ] 

-2 



a gh (q 2 



1 + 



Hi 2 ) 

1 + G(q 2 )\ 



a(q 2 ). 



(2.16) 
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An important corollary of the dynamical equations of Eq. (13. 5 j) is that L(0) = 0. In addition, 
under very general conditions, G(0) G (—1, 0). Therefore, from (12.161) one concludes that 

a gh (0) = a(0). (2.17) 

Since, finally, a sh (q 2 ) and a(g 2 ) coincide in the deep UV, where they both reproduce the 
correct perturbative behavior, the two charges can only differ appreciably in the intermediate 



region of momenta; however, since L(q 2 ) is numerically suppressed 20|], this difference is 
rather small. 



III. EFFECTIVE CHARGES FROM LATTICE 



This section contains the main results of this article, and is composed of several subsec- 
tions. After presenting a collection of lattice data, which firmly establish the IR finiteness 
of the conventional gluon propagator A(q 2 ) (in the Landau gauge) and the ghost-dressing 
function, we embark on the actual extraction of the effective charges from the lattice data, 
using the definitions and results of the previous section. The final results of all the analysis, 
carried out throughout this section, are shown on the right panel of Fig. [7J evidently, in the 
deep IR, both charges, a gh (g 2 ) and a(q 2 ), saturate at the same finite value, as predicted on 
general principles. 



A. Lattice results for the gluon propagator 

In this subsection we present some of the most relevant lattice results on the (Landau 
gauge) gluon propagator, given that it constitutes a central common ingredient of both 
effective charges. Even though in our analysis we will use only one set of lattice data (that 
of [2I), it is important to establish that various groups coincide on the qualitative behavior 
for the Green's functions in question. In Fig. [TJwe show the results for the gluon propagator 
obtained by three independent lattice groups . Although, for each group, the lattice 

spacing and the gauge group employed are different, all results have as a common feature the 
appearance of a plateau in the deep IR region, namely one of the most salient and distinctive 
predictions of the the gluon mass generation mechanism. In fact, the three set of data can 
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FIG. 1: The lattice data for the gluon propagator obtained by three independent groups 
fitted by Eq. (|3.ip . Upper left panel: Lattice data from Ref. [s]] renormalized at \i = 3.0 (AcY. 
Upper right panel: Lattice result, renormalized at p = 3.0 GeV, obtained in Ref. i|. Bottom panel: 
The SU(2) gluon propagator obtained in Ref. 



be accurately fitted in terms of a massive gluon propagator of the type 



2 i 2 

m + q 



13C A gj ( q 2 + pm 2 
In 



(3.1) 



96tt 2 v P 

where m 2 , g 2 , and p are treated as free fitting parameters. For the SU(3) lattice simulations, 
p, will be chosen to coincide with the renormalization point, while for the SU(2) case we will 
treat it as an adjustable parameter. 

Specifically, for the data presented on the upper left panel {3], we find that 
0.16 GeV 2 , g 2 = 8.79, p = 4, and p = 3 GeV. For the upper right panel, we use 



m 



m 



0.11 GeV2, g) = 9.77, p = 9.6, and p = 3 GeV; while the SU(2) lattice data of Ref. 



.a 



can be accurately adjusted using m 2 = 0.27 GeV 2 , gj = 27.68, p — 4, and p = 1.92 GeV. 

The parameter m acts as a physical mass scale, whose function is to regulate the pertur- 
bative RG logarithm; so, instead of diverging at the Landau pole, the logarithm saturates at 



a finite value. Clearly, for large values of q 2 , we recover the one-loop expression of the gluon 
propagator in the Landau gauge. Note also that for the purposes of this fit we have treated 
m as if it were a hard mass, even though an important theoretical feature of the dynamically 
generated mass is that it should be function of the momentum, vanishing in the deep UV 
in a way consistent with the operator-product expansion, displaying either logarithmic or a 
power-law running j^J [viz. Eq. fl3.6p . 

Even though it is evident from Fig. [I] that the various lattice groups appear to be in 
qualitative agreement with each other, for the actual extraction the effective charges we will 
rely on the data of Ref.^J, given that this latter group uses SU(3) simulations, and has 
available data also on the ghost propagator. 



B. The remaining lattice ingredients: numerical fits and ^-dependence 

As mentioned earlier, the different definitions of the QCD effective charges involve three 
fundamental Green's functions: the gluon propagator A(g 2 ), the ghost dressing function 
F(q 2 ), and the auxiliary function G(q 2 ). In order to verify explicitly the expected //- 
independence of r(q 2 ) and d(q 2 ), we need to have at our disposal lattice data for A(g 2 ), 
F(q 2 ), and G(q 2 ) at different renormalization points. To that end, we will exploit the 
property of multiplicative renormalizability, which allows one to connect a set of points 
renormalized at fi with the corresponding set renormalized at v, through the relations 

Using into Eq. fl3.2p the fundamental identity of Eq. (12. 15ft . whose form must be preserved 
after renormalization, we have that G{q 2 ) must satisfy 

i , 2 2N,r/2 2n I + G(q 2 , v 2 ) + L(q 2 , is 2 ) 

1 + G(q , /x ) + L(q , M ) = i - ^ V2) - . (3.3) 

Evidently, the self-consistent renormalization procedure of G(q 2 ) requires the knowledge of 
L(q 2 ). However, in Ref. 40] . the renormalization of G(q 2 ) was carried out considering L(q 2 ) 



to vanish for all values of momentum, i.e., setting L(q 2 ) = 0. To be sure, this approximation 



37] 



will not produce any appreciable difference in the deep IR region, where it was shown 
that indeed L(0) = 0; on the other hand, minor changes in the intermediate and UV regimes 
are to be expected, which, however, will be neglected for the purposes of this subsection. 
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L(q 2 ) will be eventually obtained indirectly, by substituting lattice data for F and A into 
the dynamical equation of (13.51) ; it turns out that it is indeed numerically suppressed within 
the entire range of available momenta, see Fig. 

Thus, considering for now L(q 2 ) = 0, and choosing the two different values fi = 2.5 GeV, 
and \i = 4.0 GeV, we obtain the curves for A(q 2 ), F(q 2 ), and G(q 2 ) shown in Fig. [2J 

As can be seen from this figure, the lattice sets used are restricted to momenta ranging 
roughly from 0.01 to 22 GeV 2 . Within this range, the gluon propagator can be fitted by 
Eq. (13.11) . while F(q 2 ) and G(q 2 ) are being given by 

F{q2) = TTWMr + ° 2 ' G{q2) = i + WMf " &2 ' (3 ' 4) 

with the values of the fitting parameters quoted in the caption of Fig. [2j 

Notice that in the case of the ghost dressing function, the lattice data, and correspondingly 
our fit, show no enhancement in the deep IR; instead, F saturates at the constant value 
a\ (in agreement with the large- volume lattice simulations). In addition, the Kugo-Ojima 
confinement criterion is clearly not satisfied, since G(0) deviates appreciably from the special 
value of —1. 



C. Fixing the value of g 2 (//) 



The next step is to determine the value of the renormalized coupling g(fi 2 ) that enters 
in both the definitions of the two RG-invariant quantities, r(g 2 ) and d(q 2 ). To that end, 
we resort to the three (renormalized) dynamical equations for F(q 2 ), G(q 2 ), and L(q 2 ) j^tj], 
namely 

{k-qf- 



F-\q 2 ) 

1 + G(q 2 ) 
L{q 2 ) 



Z c + g 2 C A 



k 2 q 2 



Z c + 



9 2 C A 
d- 1 



(d-2) + 



A(k)D(k + q), 



(k-q)' 
k 2 q 2 



A{k)D(k + q), 



g 2 c A 

d-lji. 



d 



(k ■ qf 
k 2 q 2 



A{k)D{k + q). 



(3.5) 



The above equations have been derived using tree-level values for the two fully dressed ver- 
tices appearing in them, namely the conventional ghost-gluon vertex and the kernel H^ u 
[viz. Eq. fl2.11j) ]. This appears to be a good approximation, given that (i) T M has been studied 
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FIG. 2: Lattice results for the gluon propagator (upper left panel), the ghost dressing func- 
tion (upper right panel) and —G(q 2 ) (bottom panel) renormalized at fi = 2.5 GeV (red curve) 
and // = 4.0GeV (black curve). Values of the fitting parameters of Eqs. ()3. 1 j) and (|3.4|) are: 



m 



m 



0.18 GeV 2 , gj = 10.64, oi = 2.64, a 2 = 0.84, b x = 0.58, and b 2 = -0.14, for // = 2.5 GeV; 
0.14 GeV 2 , = 6.95, at = 2.89, a 2 = 0.91, b% = 0.65, and b 2 = 0.047, for = 4.0 GeV. For 



both values of fj, we use p = 4, pi = 0.8, g 2 = 0.36 GeV 2 , p 2 = 0.98, and r?| = 1.45 GeV 2 . 



in lattice simulations 



43| . where it was found to deviate only mildly from its tree- level value, 



and (ii) and H^ v are connected by the identity f )2.12p . The renormalization constant Z c 
is determined by the condition F(fi 2 ) = 1. Notice that the renormalization procedure fol- 



lowed 



20( preserves the form of the crucial BRST identity Eq. f )2.15p . as required. The most 



immediate consequence of this renormalization procedure is that the value of G(fi 2 ) ^ 0; in 
fact, G(fi 2 ) = — L(/i 2 ) [see also the discussion following Eq. (13. 3p . 

The way the value of g(fJ, 2 ) is determined from Eq. (13. 5 p is the following. One substitutes 
into the integrals on the rhs the lattice data for A and D, carries out the integration 
numerically, and adjusts the value of the g 2 multiplying the integrals such that the result of 
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FIG. 3: Comparison between the ghost dressing function F(q 2 ) obtained from the ghost SDE 
(continuous lines) and the corresponding lattice data at \i = 2.5 GeV (left panel) and fi = 4.0 GeV 
(right panel). 




FIG. 4: Comparison between the G(q 2 ) function obtained from the SDE (continuous lines) and 
the corresponding lattice data at \i = 2.5 GeV (left panel) and /u = 4.0 GeV (right panel). 

the integration coincides as well as possible with the available lattice data on F and G. 

The results of this procedure for the ghost dressing function F(q 2 ) are presented in Fig. [3l 
and for G(q 2 ) in Fig. HJ in particular, we obtain the value a(/i 2 ) = 0.467 for /i = 2.5 GeV, 
and a(/i 2 ) = 0.309 for /i = 4.0 GeV. 

In order to check if the values found for a(fi 2 ) through the above procedure are compatible 
with what one would expect within the momentum subtraction scheme (MOM) that we use, 



we compare them with the corresponding four-loop perturbative calculation presented in 



44] 
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FIG. 5: The perturbative running coupling in the MOM scheme, «mom(^ 2 )i up to four-loops for 
different values of A QCD . The black squares represent the values we use for a(/u 2 ). 

The result of this comparison is shown in Fig. J2 the yellow band is obtained by varying the 



A QCD , appearing in the expression derived in 44|, in the range between 350 — 450 MeV. As 



we can see, the best adjustment for the values of a(/i 2 ) occurs for A 



QCD 



410 MeV. 



D. Final results 



From all the ingredients presented so far, one may construct the two RG-invariant quan- 
tities, r(q 2 ) and d(q 2 ) of Eqs (I2.3P and ( I2.13p . A crucial check of the self-consistency of the 
entire procedure is the numerical verification of the theoretically expected independence of 
the above quantities of the renormalization point [i. To verify this important point, r(g 2 ) and 
d(q 2 ) have been calculated using into the defining equations two different sets of inputs for 
A, F, and G, one set renormalized at /i=4.0 GeV, and another renormalized at /i=2.5 GeV. 
The values for a(/i 2 ) are precisely those obtained through the procedure of the previous 
subsection, namely a(fi 2 ) = 0.467 for /i = 2.5 GeV, and a(/i 2 ) = 0.309 for fi = 4.0 GeV. 
The results of this construction are shown on the left panel of Fig. [7J clearly, the r(g 2 ) and 
d(q 2 ) obtained from each set of data are practically on top of each other, thus numerically 
confirming the theoretical expectations. One can also see that the two quantities behave as 
expected, differing only in the intermediate region of momenta (20 - 600 MeV). 

Next, using as ingredients the lattice data for A(q 2 ) and F(q 2 ) presented in Fig. [2] and 
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FIG. 6: L(q 2 ) determined from the corresponding SDE (|3.5p . using the solutions for A(g 2 ) and 
F(q 2 ) presented in the Fig. [2] respectively, at the same renormalization point. 

the values of a(/i 2 ) quoted in Fig. one can compute, for the sake of completeness, the 
auxiliary function L(q 2 ) from the SDE (13. 5p . The results for L(q 2 ) are presented in Fig. |6]for 
both values of fi 2 . From Fig. |6l it is easy to check the three properties of L(q 2 ) mentioned 
before: (i) indeed L(q 2 ) is numerically rather small over the full range of momenta, (ii) it 
vanishes in the deep IR, and (iii) its maximum occurs in the intermediate momenta region 
(around 500 MeV). 

At this point, the non-perturbative running charges, a gh (q 2 ), and ct{q 2 ), defined in 
Eqs. (I2.4p and (I2.9p . respectively, may be extracted by multiplying the results obtained 
for the corresponding RG-invariant quantities by the factor [q 2 + m 2 (q 2 )]. To do that, how- 
ever, one must assume a functional form for the running mass m 2 (q 2 ). We will use a mass 
that decreases in the UV as power-law running (see, e.g., [42|, |45| ) 

m 2 (g 2 ) = KK<f + m l) • (3-6) 

The running mass of (I3.6P has a finite value at q 2 — > 0, i.e. m 2 (0) = m 2 ,, with a power-law 
decrease in the deep UV. For m we choose some representative values consistent with the 
phenomenological studies, namely m = (500 — 600) MeV jsl, l^fj (47). 

The effective charges obtained following the above steps are shown in the right panel of 
Fig. [71 Evidently, both charges exhibit the correct (UV) perturbative behavior, and freeze 
at the same finite IR values corresponding to a gh (0) = a(0) = 4.45 (m = 500 MeV) and 
Q! gh(0) = ct(0) = 6.40 (m = 600 MeV). The difference between the two couplings is only in 
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FIG. 7: Left panel: Comparison between the two RG-invariant products d(q 2 ) (solid line) and 
r(q 2 ) (dashed line); notice that there are two overlapping curves at different fi for each product. 
Right panel: Comparison between the QCD effective charge extracted from lattice data: a(q 2 ) (red 
line with circles) and a gh (g 2 ) (black line with squares) for two different masses: too = 500 MeV 
(dashed) and to = 600 MeV (solid). 

the intermediate momenta region, and it is entirely due to the L(q 2 ) function; in this region 
a(q 2 ) is always bigger than a gh (q 2 ). 

IV. RECONCILING LATTICE WITH PHENOMENOLOGY 

The effective charges we have obtained from the lattice (within the MOM renormalization 
scheme that we use), reach values in the deep IR that are almost an order of magnitude 
higher than those obtained from a large number of phenomenological studies. In particular, 
while the charge obtained from the lattice ranges between 4.5 - 6.5 [depending on the value 
of m(0)], the systematic fitting of numerous processes suggests values for the effective charge 
in the range 0.7 ± 0.3 (for a similar range of gluon masses). 

If one were to take both lattice and phenomenological results at face value, one should 
attempt to determine the reason for this sizable discrepancy. In this section we will address 
this issue in the context of a toy calculation, and we will argue that the observed discrepancy 
may be traced back to the difference in the gauge used when extracting the lattice results 
(the Landau gauge of the BFM, £ Q = 0) and that assumed in the phenomenological studies 
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(the Feynman gauge of the BFM, £ Q = 1). Even though we cannot reach firm conclusions, 
our calculation seems to indicate that the difference in the gauges may indeed reconcile 
lattice with phenomenology. 

The crux of the matter is that the effective charge entering into physical processes is 
neither a gh (q 2 ) nor a(q 2 ), but rather the genuine PT effective charge js-8], to be denoted by 
a PT (g 2 ). This charge is defined exactly as a(q 2 ) in ( 12. 9p . but with the crucial difference that 
the propagator used to form the d(q 2 ) is the PT gluon propagator, i.e., the BFM propagator 
calculated in the Feynman gauge. The Feynman gauge of the BFM is privileged, in the 
sense that it is selected dynamically when the gluon self-energy is embedded into a phys- 
ical observable (such as an on-shell test-amplitude). Specifically, this gauge captures the 
net propagator-like subamplitude emerging after QED-like properties have been replicated 
inside the test-amplitude, by means of the PT procedure. Therefore, any gauge-related 
exchanges between the Green's functions put together to form observables, are eliminated 
in this particular gauge. Instead, the gluon propagator in the Landau gauge, for example, 
contains still residual unphysical contributions, which, when introduced into a physical am- 
plitude, will cancel against similar terms from vertex and box diagrams (see third item in 
3). 

For asymptotically large momenta the numerical difference between the charges defined in 
either gauge is controlled by the constant C^ Q , given in (12. 6p . Evidently, in the UV this dif- 
ference is subleading, and cannot give rise to any appreciable difference. Non-perturbatively, 
however, the difference between the two charges may be sizable. If, for example, we subscribe 
to the notion of dynamical mass generation, a difference in the gauge may lead to a vastly 
different IR behavior. In order to gain a quantitative understanding of how the difference in 
the gauge used may cause a significant disparity in the infrared values of the corresponding 
effective charges, we consider a model where the gauge bosons are endowed with a mass at 
tree-level. This will allow us to calculate, at one-loop level, the deviation between the two 
propagators, and the discrepancy that it induces to a(0) and a PT (0). The model in question 
is simply the electroweak sector of the Standard Model, with the electric charge set to zero, 
or, equivalently, with sin0 w = 0, where 8 W is the electroweak mixing (Weinberg) angle. In 
this limit the three gauge bosons (two Ws and one Z) are degenerate. 

At one-loop, the SU(2) gluon self-energy, or equivalently, the Z-boson self-energy, may 
be obtained from the results of 4f|, for any value of £ Q . Specifically, one has (in Minkowski 
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space) 



- [4 (5 - 2d) m A - 8 (2 - d) m 2 g 2 + m 2 g 2 + g 4 ] ^(g 2 ; m 2 , m 2 ) 

+ 2^-^ [(1 - £ Q ) 2 m 4 - 2 (3 - 2d + £ Q ) m 2 g 2 + g 4 ] £ (g 2 ; m 2 , £ Q m 2 ) 

- (g 2 + 5m 2 ) (g 2 -4£ Q m 2 )£ (g 2 ;£ Q m 2 ,£X), (4.1) 

where IlF(g 2 ) = II(g =1 )(g 2 ) is the conventional gluon self-energy in the Feynman gauge, m 
denotes the effective gauge boson mass, and A and B are given by 

A (m 2 ) = 16vr 2 ' ' 



k 2 



,2 



^(g 2 ;/?! 2 ^ 2 ) = 16tt 2 / — 1 ^ . (4.2) 

A {k 2 ~ mf)[(k + q) 2 - mi] 

Setting £ Q = 1 in the above formula we recover the standard PT result for the one-loop 
self-energy of the Z-boson {49| , to be denoted by flF(g 2 ), namely 

n F (g 2 ) = n F (g 2 ) - £^(g 2 - m 2 )B (q 2 ;m 2 ,m 2 ) . (4.3) 

Let us now take the difference R(q 2 ) between U^ Q (q 2 ) calculated in the Landau and Feynman 
gauges (£ Q = 0, and £ Q = 1, respectively); denoting the former by n L (g 2 ), one has in the 
limit g 2 — > and d = 4, 



R(o) = n L (o)-n F (o) 

( m 4 d 



y 2. ^2 ^2 



S (g 2 ; m 2 , 0) + 3 [5 (g 2 ; m 2 , 0) - B (q % ; m~, rrr 



■ (4-4) 

q 2 =0 



(An) 2 { 6 dq 2 

We next extend the one-loop expression given in (14. 4 p to the non-perturbative regime, by 
introducing the following approximations: (i) we replace the (tree-level) massive propagators 
appearing in the function B by their fully dressed counterpart A (in Landau gauge), and 
(ii) the "hard" mass m by its running counterpart. Then we find (in Euclidean space) 



d_ 
dq 



2 5 (g 2 ;m 2 ,0) 



B (q 2 ;m 2 } 0)-B (q 2 ;m 2 } m 2 )\ 



-> — 

7T- 



1 d f A(k 2 ) 
7T 2 dq 2 J k (k + q)< 

m 2 (k 2 )A 2 (k 2 ) 
k 2 



:A(0), 



(4.5) 
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A(0) [GeV" 2 ] 


/i [GeV] 




m [MeV] 


q(0) 


R(0) 


«pt(0) 


5.51 


2.5 


0.467 


600 


6.40 


0.31 


0.47 


7.00 


4.0 


0.309 


600 


6.40 


0.33 


0.30 


5.51 


2.5 


0.467 


500 


4.45 


0.15 


0.62 


7.00 


4.0 


0.309 


500 


4.45 


0.16 


0.40 



TABLE I: The gauge-invariant and universal IR fixed point a PT (0) obtained from the Landau 
gauge a(0) one after applying Eq. (|4.8p . 



Thus, in the SU(3) case one obtains the final result (with y = k 2 



R(0) 



3 a(fj, 2 )ml 
2 4V 



2 r 



mo 
12 



A(0)+3 dym 2 (y)A 2 (y) 



(4.6) 



where the multiplicative factor of 3/2 corresponds to the ratio of the Casimir eigenvalues 
for the adjoint representations of the gauge groups SU(3) and SU(2). 
Since in Euclidean space R(q 2 ) changes sign, we obtain 1 

1 1 



q 2 + U F (q 2 ) 
A(q 2 ) 



(q 2 + U L (q 2 )) 1 + 



R{q 2 ) 



1 + R{q 2 )A{q 2 ) 
arriving at the following relation for the two couplings, 

a(0) 



(4.7) 



a PT (0) 



1 + J R(0)A(0) 

In order to get an approximate estimate for ct PT (0) we need to determine the value of 
i?(0) from (14.61) . To that end, we use the lattice data for the A(y) appearing on the rhs, and 
a mass m 2 (y) that displays power-law running, given by (I3.6p . The results of this procedure 
are summarized in Table HJ clearly, the values obtained for a PT (0) are indeed much closer to 
the expectations based on phenomenological studies. 

We emphasize that Eq. (14. 61) constitutes only a simplified estimate of the complete answer, 
and our results are suggestive at best. Note in particular that, as is evident from Table fl], 



1 To go to Euclidean space, we set q 2 = — g|, with q 2 , > the positive square of a Euclidean four-vector, 
define the Euclidean propagator as A E (q 2 ) = — A(— q 2 ,), and the integration measure as J k — i j k ^. To 
avoid notational clutter we always suppress the subscript "E" . 
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Eq. (I4.8p leads to the introduction of a spurious dependence on the renormalization scale \i 
for the ostensibly RG- invariant quantity a PT (q 2 ). 

V. CONCLUSIONS AND OUTLOOK 

In this article we have shown how to extract effective QCD charges from the available 
(quenched) lattice data for some of the fundamental Green's functions of QCD. We use two 
different definitions of the effective charge, whose construction follows a similar procedure, 
relying on the construction of RG- invariant quantities out of the judicious combination of the 
various field-theoretic ingredients. The effective charges obtained display the characteristic 
feature of freezing at a common finite (non-vanishing) value in the deep IR, as expected 
from a variety of theoretical and phenomenological considerations. 

In addition, we have offered a plausible explanation for the observed discrepancy in the 
freezing values of the effective charges obtained from the lattice and those derived from the 
fitting of various QCD processes, sensitive to non-perturbative physics. Our claim is that the 
underlying reason for the discrepancy is the difference in the gauges (Landau vs Feynman) 
used in the two approaches. We have studied this issue in the context of a toy model, which 
seems to corroborate this assertion. 

It is clearly highly desirable to have available lattice results for the gluon and ghost 
propagators in gauges other than the Landau. In fact, a new gauge-fixing algorithm that 
may all^ one to carry out .attics stations in genera. R t gauges has been recent.y pro- 
posed [50!]. In addition, it is of considerable theoretical importance to obtain lattice results 



in the Feynman gauge of the BFM [5l|], where, by virtue of the PT, quantities such as the 
gluon propagator acquire a gauge-invariant and universal status. Lattice results in this class 
of gauges would allow not only a direct determination of the phenomenologically relevant 



coupling a PT (0), but will furnish a stringent test of the SDE predictions for the gluon 41] 
and ghost propagators |52[ ] . 
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